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ABSTRACT 

Aims. The objective of this work is to obtain an extinction-corrected distribution of optical surface brightness and colour indices of 
the large nearby galaxy M 3 1 using homogeneous observational data and a model for intrinsic extinction. 

Methods. We process the Sloan Digital Sky Survey (SDSS) images in ugriz passbands and construct corresponding mosaic images, 
taking special care of subtracting the varying sky background. We apply the galactic model developed in Tempel et al. (2010) and 
far-infrared imaging to correct the photometry for intrinsic dust effects. 

Results. We obtain observed and dust-corrected distributions of the surface brightness of M 3 1 and a map of line-of-sight extinctions 
inside the galaxy. Our extinction model suggests that either M31 is intrinsically non-symmetric along the minor axis or the dust 
properties differ from those of the Milky Way. Assuming the latter case, we present the surface brightness distributions and integral 
photometry for the Sloan filters as well as the standard UBVRI system. We find the following intrinsic integral colour indices for M 3 1 : 
(U - B) = 0.35; (B - V) = 0.86; (V - R) = 0.63; (R - I) = 0.53; the total intrinsic absorption-corrected luminosities of M31 in 
the B and the V filters are 4.10 and 3.24 mag, respectively. 



Key words, galaxies: individual: Andromeda (M31) - 
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y- 1. Introduction 

The primary information about the structure of galaxies and 
,— I ' their stellar populations comes from the distribution of sur- 
| face brightnesses and colour indices. Thus, it is only natu- 
ral that the photometry of the largest nearby galaxy M 3 1 has 
bee n studied thoroughly for almost a century (for references , 
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see Ivan den BergUl99UTTenies et al.lll994tlTempel et"alll2010l) . 
More recently, satellite observations in the far-infrared have 
become available and attempts have been made to consider 
the i mpact of dust extinction on the observable properties of 
M 31 (IXu &Heloul 19961: lHaas et al.|[T998tlMontalto et al.l2009t 
Tem pel et al.ll2.Q10» . While the general luminosity distributions 
of different studies are in good agreement, the distribution of 
colours is not so well settled. Small systematic deviations be- 
tween different datasets lead to considerable uncertainties of 
colour indices - too big for studying in detail the stellar pop- 
ulations or the star formation history of the galaxy with the help 
of chemical evolution models. 

Recently, w ithin the Sloan Digital Sky Survey (SDSS; 
lYork et al.l l2000l) . a contiguous strip covering the entirety of 
M 3 1 through ugriz filters has been observed. The combination 
of five filters on a single, well-calibrated telescope and roughly 
uniform atmospheric conditions during the observations provide 
a uniquely homogeneous dataset for a thorough analysis of the 
extensive galaxy. In this paper, we use the SDSS observations 
to study the detailed luminosity and colour distribution of M 3 1 . 
We use also far-infrared imaging by the Spitzer Space Telescope 
and Infrared Astronomical Satellite to correct the derived pho- 
tometry from dust extinction. 
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Fig. 1. Upper panel (a): SDSS imaging of M 31 in g filter, cor- 
rected for background variations. Middle panel (b): the same 
as upper panel, with intrinsic extinction effects removed (see 
Fig- Eh for extinction map). Lower panel (c): the same as middle 
panel, showing also the masked regions (white circles). 



Thro ughout the paper, we assu me the distance of M 3 1 to be 
785 kpc dMcConnachie et a"Dl2005l) . corresponding to the scale 
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Fig. 2. Elliptically averaged observed surface brightness profiles from the SDSS photometry, converted to the UBVR I system (black 
solid lines) together with uncertainties (grey regions). Conversion to the UBVRI system is done according to iBlanton & Rowels! 
(l2007h . For comparison, some earlier measurements are also shown (points with error bars; for references, see Sect.[H> . The earlier 
measurements in UBVR are from elliptically averaged data, while / and outermost parts of V measurements by llrwin et al.l (l2005h 
have been made along the minor axis only. Left panel: major axis; right panel: minor axis. For clarity, the BVRI profiles have been 
shifted by the indicated values. 



Table 1. The applied solar luminosities (in AB magnitudes for 
ugriz filters and in Vega magnitudes for UBVRI filters) and 
Galactic extinctions (A) for each filter. 
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1 ' = 228 pc. The inclination angle of the galaxy has been take n 
1T.5 (IWalterbos & Kennicutdll988tlde Vaucoule urs et al.lll99lh 
Accor ding to IWalterbos & Kennicuttl (11987 ). Ferguson et alj 
(12001 . and our analysis of the SDSS and Spitzer images, we 
set 38° 1 as the major axis position angle. All luminosities 
and colour indices have been corrected f rom e xtinction in the 
Milky Way according to Schleg el et al.l ([1998) and are given 
in AB-magnitudes for the ugriz filters and in Vega magnitudes 
for the UBVRI filters, as usual. For the Sloan filters, the extinc- 



tion is derived from the Schl egel et al.l (1 19981) estimates and the 
Galactic extinction law by linear interpolatio n. The absolute so- 
lar lum inosity for each filter was taken from Bl anton & Roweisl 
(2007). The applied solar luminosities and the Galactic extinc- 
tion for each filter are presented in Table [T] 



2. SDSS data reduction 

M 3 1 was observed with the Sloan telescope on October 6, 2002, 
using ugriz filters. We retrieved the pipeline-corrected frames 
(fpC-files) from the SDSS Data Archive Server (DAS). The 
frames (2048x1489 pixels, 0.396 arcsecpx 1 ) originate from 
two observing runs obtained with six columns of CCDs, thus 
giving a total of 12 scanlines parallel to the apparent major axis 
of the galaxy. The scanlines have a width of 13.52 arcmin, over- 
lapping by about 55 arcsec. We used 56 frames from each scan- 
line (field numbers 49-104), in order to construct mosaic images 
in ugriz passbands. 



Individual frames were aligned by applying astrometric 
transformations given in frame headers. In the final mosaics, 
10x10 binning was implemented, resulting in image dimensions 
of 2300x7400 pixels. For the central region of M 31, additional 
images with higher sampling were created using 2x2 binning. 

Sky background was estimated by analysing the SDSS 
frames in several steps. We used the sky.pro routine from the 
IDL Astronomy User's LibrarjQ which gives an estimated mode 
of the sky-background values. The total sky contribution in each 
pixel, sfcytotai) was assumed to consist of three components: 



S^ytotal — S^yscanline + Subtemporal + sky gradient- 



(1) 



The base sky value, sky scw \i ne , was determined at both ends of 
each scanline, at a sufficient distance from M 3 1 . Superimposed 
on it, variations along and perpendicular to the scanlines were 
found, probably caused by different effects. Most likely, the for- 
mer is the result of temporal sky brightness variations during 
the observations and is thus designated here as sky tempora \. It was 
determined on the basis of the two outermost scanlines, sep- 
arately for both observing runs, with a moving window with 
the width of a full frame and the height of 1000 pixels, yield- 
ing a smooth curve along the scanline. Sky variations perpen- 
dicular to the scanlines (and perpendicular to the major axis 
of M31), sky gradient, were probably caused by stray light from 
the galaxy itself. Its contribution was determined by matching 
the sky-subtracted signal in overlapping regions of side-by-side 
scanlines. The values of sAry gra dient were then linearly interpolated 
between the start and the end of each pixel row of each scanline. 

After subtracting the soft bias of 1000 ADU (added in the 
SDSS pipeline) and the sky background, the pixel values were 
corrected for atmospheric extinction and were converted to 
the flux and magnitude scale, following the SDSS photometric 
flux calibration^ In order to eliminate the contribution of fore- 
ground stars and the galaxies M 32 and NGC 205, circular masks 
were applied on top of these objects (see Fig. QJ). SExtractor 



http://idlastro.gsfc. nasa. gov/ 

http://www.sdss.Org/dr7/algorithms/fluxcal.html#counts2mag 
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Fig. 3. Surface brightness profiles along the minor axis before 
and after intrinsic dust corrections. Red dotted line shows the 
profile along the nearer side (NW direction); blue dashed line 
shows the profile along the further side (SE direction). The cor- 
responding elliptically averaged surface brightness profiles are 
also shown (black solid lines) together with error estimates (grey 
region). For clarity, the dust-corrected profiles have been shifted 
up by two units. 
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Fig. 4. Panel (a): map of the total optical depth r max for the Sloan 
g filter; panel (b): extinction map for the g filter, calculated ac- 
cording to the best-fit model; panel (c): the elliptically averaged 
profile of the optical depth map for the g filter; grey region indi- 
cates the scatter of the optical depth along each ellipse. 



dBertin & Arnoutsl 1996b was used to detect and mask faint back- 
ground galaxies in the outer regions. 

For the uncertainty estimation, the following noise sources 
were taken into account: the random photon noise, the contri- 
bution of dark current and readout noise, the uncertainty of the 
sky level estimate (0.5 ADU), and the error of matching bright- 
ness levels in the overlapping regions (0.5 ADU). More details 
about the conducted image processing and sky removal will be 
provided in a separate paper (in preparation). 

The final image of M 31 for the g filter is shown in Fig.QJ. 
The applied masks are shown in Fig.QJ. 

3. Structural model and extinction corrections 

To correct the optical images for dust extinction, we have to con- 
sider the spatial distribution of stars and dust inside the galaxy. 
For this purpose, we have relied on a general three-dimensional 
galaxy model and extinction calculation meth ods, briefly re - 
viewed below. More details can be found in Tem pel et al.l d2010). 

The model galaxy is given as a superposition of its individ- 
ual stellar components and a dust disc. Each stellar component 
is approximated by an ellipsoid of rotational symmetry with 
a constant axial ratio q; its spatial density distribution follows 
Einasto's law 



1(a) = Iq exp 



a 
kao 



i/N 



(2) 



where Iq = hL/fAnqai) is the central density and L is the compo- 
nent luminosity; a = ^r 2 + z 2 Iq 2 , where r and z are two cylin- 
drical coordinates; an is the harmonic mean radius. The coeffi- 
cients h and k are normalising parameters, depending on N. 

In the case of the young disc the spatial density distribution 
is assumed to have a toroidal form, approximated as a super- 
position of a positive and a negative density component, both 
following Eq. [2] 



Density distributions of all stellar components are projected 
along the line of sight and their sum yields the surface brightness 
distribution of the model 
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(3) 



(4) 



where l(r,z) denotes the spatial luminosity density (Eq. [2J of the 
galaxy component, L(X, Y) is the corresponding surface bright- 
ness distribution, X and Y are coordinates in the plane of the 
sky, T(r, z) is the optical depth of dust along the line of sight for 
the given location, i is the angle between the sky plane and the 
galaxy plane, and the summation index k designates each visible 
component. 

A dust disc is added to the model galaxy, with its vertical 
density distribution following Eq. [2] Along each line of sight, 
the optical depth r(r, z) is assumed to be proportional to the dust 
column density. The optical depth is taken zero between the ob- 
server and the dust disc, and equal to r max (X, Y) behind the dust 
disc. T max is the total optical depth along a given line of sight and 
is thus a function of (X, Y). Inside the dust disc, r(r, z) i s between 
and r max ; see fig. 1 and eq. 6 in Tempel et al.l (12010b . 

Assuming that dust column density is proportional to the 
far-infrared flux, the map of r max (X, Y) can be derived from far- 
infrared imaging: 



'max,/ 

(X, Y) = c f j,j,F A ,(X, Y) 



A^B(A', T ) 
(A'fB(A, T(X,Y))' 



(5) 



where F^(X, Y) is the far-infrared flux map at a reference wave- 
length A'; T is a reference temperature; T(X, Y) is a map of the 
dust temperatures; B(A, T) is the black-body function; B is the 
dust emissivity index; and c is an empirical calibration constant, 
corresponding to the reference wavelength and temperature, and 
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to the filter / of the optical observations (c oc A^IE{B - V) = 
RvAa/A v ). 

We used the Spitzer Space Telescope observations at 24, 70 
and 160 //m and the Infrared Astronomical Satellite detections at 
100 /mi to derive maps for the far-infrared flux and dust temper- 
ature. The point-spread-functions and sampling of the Spitzer 
observations were reduced to that of 160 /im imaging. A dust 
model, comprising a colder and a warmer component was used 
to approximate the measured far-infrared spectral energy distri- 
bution along each line of sight with modified Planck functions; 
extinction was ascribed to the colder component only. More de- 
ta ils about this step an d the related error estimates can be found 
in lTempel etai1(l2010l) . 

F or applying Eq [5] we c hose A' = 100 /mi a nd T _ 18.2K 
as in Schlegel et~atl (U998I and according to IVlahakis et al.l 
(2005) we take the dust emissivity index B = 2.0. The best value 
for the calibration constant c was found during the fitting pro- 
cess. Because of the specific geometry and viewing angle of 
M31, the surface brightness distribution along the minor axis 
is sensitive to the optical depth. Extinction is higher along the 
closer side of the minor axis because most of the bulge stars 
are behind the dust disc, whereas for the f arther side, the bulg e 
is in front of the dust disc (see fig. 1 in iTempel et al.l 12010). 
However, assuming the galaxy to be symmetric, dust-free lumi- 
nosity distributions along the minor axis in the two opposite di- 
rections should overlap, providing a recipe for the calibration of 
the constant c. During the iterative fitting process (see Sect. [4] 
and Fig. |3J, the best overlap was reached with c = 0.022 for the 
g passband. The Galactic extinction law was used to calculate 
optical depths for observations with other filters. 
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Fig. 5. Top panel shows the final model (solid lines) and ob- 
served (grey regions) profiles along the major axis in ugriz fil- 
ters. For clarity, the r, i, z profiles have been shifted by the in- 
dicated values. The errors for the model are mainly due to the 
errors of the observations. Bottom panel shows the observed ax- 
ial ratio q. The errors from the ellipse fitting for axial ratio q are 
in order of 0.01. 



4. Results: model fitting and surface brightness 
profiles 

On the basis of the SDSS observations, we have derived 
the surface brightness distributions of M31 in ugriz filters. 
Figure [2] presents the elliptically averaged observational sur- 
face brightness profiles with error estimates in UBVRI fil- 
ters. The elliptically averaged profiles were generated us- 
ing the IRAF/STSDAS task ellipse. The conversions from 
the Sloan filters to the UBV R I sys tem were conducted ac- 
cording to iBlanton & Roweisl (120071) . For comparison, some 
earlier measurements from the literature are also s hown 
(data points with error bars): Ide VaucouleursI (1 19581) (B\ 
iHoessel & Melnickl dl980t) (UBV), iHodge & Kennicuttl dl982h 
(B). iKenj (Il987l) fflUW alterbos & Kennicutj dl988l) ( UBVR) 
Pritc het& van den Berghl (11994 (V). and llrwin et all (120051) 
(W). The UBVR profiles are elliptically averaged profiles, the / 
profile is measured along the minor axis of M 3 1 . 

However, for applications in stellar population analyses, sur- 
face brightness and colour distributions corrected for intrinsic 
dust extinction are more desirable. The amount of extinction de- 
pends on the properties of the dust disc and on the general three- 
dimensional luminosity distribution of a galaxy. The correspond- 
ing model is briefly described in Sect. [3] 

To fit the model to the observations, we have assumed the 
galaxy M 3 1 to consist of a bulge, an old and a young disc, an 
outer halo, and a central nucleus. Mo del construction deta ils are 
similar to the modelling presented in Tempel et aQ (l2010l) . 

For model fitting we used the Sloan ugriz profiles. Model fit- 
ting was started with crude initial model parameters. In the first 
step of the fitting process, a stellar model was fitted to the ellip- 
tically averaged surface brightness profiles (using ellipse), ne- 
glecting the dust disc and thus any intrinsic extinction. Next, the 



dust disc was added and the extinction map was calculated ac- 
cording to the optical depth map and the luminosity density dis- 
tribution model, derived in the first step. Now the extinction map 
was applied to the full 2-dimensional optical imaging, yielding 
a dust-free image of the galaxy. The calibration constant c in 
Eq. [5] was determined, providing the best overlap of the dust- 
free surface brightness profiles along both sides of the minor 
axis. In Fig.[3]the profiles along both sides of the minor axis be- 
fore and after dust corrections are shown. Finally, the structural 
model was fitted to the extinction-corrected imaging, allowing 
us to constrain the model parameters. All these steps were re- 
peated iteratively: according to the new model, a new extinction 
map was calculated and applied to the initial imaging; calibra- 
tion constant c was adjusted and a new extinction-corrected im- 
age was calculated, allowing us to further constrain the model 
parameters. Four iterations were needed to reach convergence. 
The final model is presented in AppendixlAl 

Figure|4]shows the total optical depth map (panel a) and the 
corresponding extinction map (panel b) of M 3 1 for the Sloan 
g filter. Extinction is higher along the closer side of the minor 
axis due to the geometry and viewing angle of the galaxy. This 
asymmetry was used to determine the calibration constant c in 
Eq. |5] during the fitting process as described above. Figure HJ; 
presents the elliptically averaged profile of the total optical depth 
for the g filter; the grey region illustrates deviations of the optical 
depth values along each ellipse. In Fig.|4j;, the signature of three 
dust rings can be seen. 

Figure [T] displays the background-subtracted ^-filter mosaic 
images of M 3 1 before (panel a) and after (panel b) dust extinc- 
tion corrections. The profiles presented in Fig. [3] are measured 
from these images. 
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Table 2. Visible and intrinsic colour indices of M31, corrected 
for extinction in the Milky Way. 
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0.43 


0.86 






4.36 



Notes: {fl) estimated errors are 0.09 mag for u-g and 0.06 mag 
for other colour indices. esti mated errors are on the orde r of 
O.lOmag. Colour indices given by de Vaucouleurs et al. (1991) are 
visible colour indices. 




Radius (major axis) (kpc) 

Fig. 6. Distribution of the elliptically averaged observed (grey 
regions) and modelled (solid lines) colour indices for the Sloan 
filters along the major semi axis. The scatter of the grey region 
indicates the uncertainty of observations and ellipse fitting. The 
colour indices are corrected for the absorption in the Milky Way 
and for the intrinsic absorption of M 3 1 . 



The derived intrinsic ugriz surface brightness profiles along 
the major axis of the best model are shown in Fig. [5] The lower 
panel of Fig. [5] shows the axial ratio q of elliptically averaged 
profiles of M 3 1 . 

In Fig. |6l the distributions of the colour indices of the 
model (solid lines) are compared to the observations (grey re- 
gions) for the Sloan filters. In Fig. [Vj the same is shown for the 
standard UBVRI filters. In the latter case, the data points with 
error bars are observational data from the literature as com- 
piled in iTempel et al.l ([20 10). The colours of M31 are given 
in Table [2] The colours are corrected for the extinction in 
the Milky Way. For compa rison, the colour indices derived by 
Ide Vau couleur s et alj d!991l) are given. 

The derived value for the con stant c _ 0.022 can b e com- 
pared to the constant p in eq. 22 of lSchlegel et al.l d!998l) for the 
Milky Way. In case of the Galactic extinction law, p = 0.016 
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Fig. 7. Distribution of the elliptically averaged observed (grey 
regions) and modelled (solid lines) colour indices for the UBVRI 
filters along the major (four top panels) and minor (bottom 
panel) semi axis. The scatter of the grey region indicates the 
uncertainty of observations, ellipse fitting and conversion to the 
UBVRI system. Points with error -bars are observational data as 
compiled in Temp el et alj (|2010|) . shown for comparison; the / 
data were not available for the major axis. 



would transform to c - ln(10)/2.5 x A g /(A B - A v ) x p = 0.056 
for the Sloan g filter. Using this value in the case of M31, the 
restored intrinsic brightness along the minor axis on the side 
closer to us would be higher than on the farther side, because 
the closer side is more affected by dust absorption, as explained 
in the last paragraph of Sect. [3] Therefore, our model suggests 
that either M 3 1 is intrinsically non-symmetric, with the closer 
side being brighter than the farther side at respective distances 
from the centre, or the dust properties of M 3 1 differ from those 
of the Milky Way, leading to a different value of c. The former 
scenario is not sup ported by the Spitzer Space Telescope obser- 
vations at 3.6/im dBarmbv et al.l 12006). where the stellar light 
should be least affected by absorption. The latter scenario re- 
mains to be confirmed or invalidated with future observations of 
M 3 1 with the Herschel Space Observatory, providing higher an- 
gular resolution in far-infrared than the Spitzer Space Telescope, 
and a coverage of longer wavelengths. 

The derived luminosity profiles and the resulting distribu- 
tions of colour indices are based on a homogeneous dataset. 
Therefore, possible systematic errors in the colour profiles are 
smaller than in profiles derived by averaging data from different 
instruments. The improved spatial coverage and lower uncertain- 
ties provided by the SDSS data, demonstrated by Fig. [7] allow 
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us to obtain tighter constraints for the star formation history and 
the overall origin of M3 1 . 
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Appendix A: Structural model of M 31 

The derived extinction distribution assumes a certain three- 
dimensional light distribution model, consistent with the ob- 
served surface brightness distributions along the major and mi- 
nor axis of the galaxy. 

The final model parameters are given in Table IA.1I These 
parameters correspond to the dust-free galaxy model, i.e. they 
represent the intrinsic stellar emission of the galaxy. The model 
given in Table lATI is not unique: the same total surface bright- 
ness profiles can be achieved using slightly different component 
parameters. To determine the real components of the galaxy, 



Table A.l. Calculated parameters of the photometric model. 
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Nucleus 


0.0008 


0.0018 
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Notes: Structural parameters for the young disc are given for the 
positive component only. <b) The luminosity of the young disc is its 
true luminosity (L = L+ + L_, where L_ = -0.5L+). The luminosi- 
ties are corrected for the intrinsic extinction and the extinction in 
the Milky Way, and are given in the units of 10 10 Lq. 
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Fig. A.l. Extinction-corrected surface brightness distribution of 
M 3 1 along the major (upper panel) and minor (lower panel) axis 
in Sloan g filter. Individual stellar components are given with 
dashed lines. 

metallicity determinations and galaxy dynamics should be taken 
into account. For the present study, however, realistic deviations 
from the given model would cause negligible effects. 

In Fig. IA. li the contributions by each galactic component to 
the surface brightness distribution in the g filter are compared 
to the profile measured from the corresponding dust-corrected 
image. 



